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Abstract 

Possibility of asymmetric square convection is investigated numerically using a 
few mode Lorenz-like model for thermal convection in Boussinesq fluids confined 
between stress free and conducting flat boundaries. For relatively large value of 
Rayleigh number, the stationary rolls become unstable and asymmetric squares ap- 
pear as standing waves at the onset of secondary instability. Asymmetric squares, 
two dimensional rolls and again asymmetric squares with their corners shifted by 
half a wavelength form a stable limit cycle. 

PACS numbers: 47.20.Ky, 47.27.-i 



Two dimensional stationary roll-patterns are known to be the only stable solutions 
at the onset of thermal convection in a thin layer of Boussinesq fluid confined between 
conducting boundaries [1] except in the case of fluids with vanishingly small Prandtl num- 
ber [2]. The stationary convection in the form symmetric square cells are unstable [1] 
except in fluids of very high Prandtl number a at very high values of Rayleigh number 
R [3]. Assenheimer and Steinberg[4] observed an interesting possibility of hexagonal con- 
vective cells with both up- and down- flow in the center of the cell at R approximately 
twice the value of critical Rayleigh number R^.. The possibility of dual types of hexag- 
onal convective cells was also established by Clever and Busse [5]. Recently, Busse and 
Clever[6] predicted a new possibility of stationary convection in the form of asymmetric 
square cells with both up- and down- flow in the center for rigid and thermally conducting 
horizontal boundaries. 
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In this work we construct and study a Lorenz-like model[7] that describes convection 
patterns in the form of rolls and squares, both symmetric as well as asymmetric ones, in 
a thin layer of Boussinesq fluid of moderate and high Prandtl numbers {a > 2) at high 
values of R. The nonlinear superposition of mutually perpendicular sets of rolls of the same 
wave number would be called asymmetric squares when the structure (e.g., shadowgraph 
picture) does not possess four-fold symmetry. It happens when the intensities of the 
two set of rolls are different. The model allows us to study competetion between rolls 
and asymmetric squares. In the following we shall derive the model from hydrodynamic 
equations and investigate the model numerically for possible stable solutions. We then 
present the results and discuss them. 

We consider an infinite layer of Boussinesq fluid of kinematic viscosity v, ther- 
mal diffusitivity k and thickness d confined between two perfectly conducting horizontal 
boundaries and heated underneath. Using the length scale d, the time scale ^, and the 
temperature scale as the temperature difference AT between lower and upper boundaries, 
the non-dimensional form of hydrodynamic equations in Boussinesq approximation read 

9tV + (v.V)v = -Vp + a{9es + V\) (1) 
V.v = (2) 
dte + Y.ve = V^e + Rv.es (3) 

where v = (fi,f2,f3) are the velocity fields, 6 the deviation from the conductive tem- 
perature profile and p the deviation from static pressure of conductive state due to con- 
vective instability. Prandtl number a and Rayleigh number R are defined as cr = ^ and 
R^^i^ns^, where a is the coefficient of thermal expansion of the fluid, g the accelera- 
tion due to gravity. The unit vector 63 is directed vertically upward, which is assumed to 
be the positive direction of x^-axis. The boundary conditions at the stress-free conducting 
flat surfaces imply 9 — V3 = d^^v^ at X3—O, 1. Taking curl twice of the momentum 
equation (Eq. 1) and using the continuity condition (Eq. 2), the equation for vertical 
velocity reads 

dtV^s = aVSs + aVH^e - eg.V x [(a;.V)v - (v.V)w] (4) 
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where a; = Vxv is the vorticity, and V — du + ^22 is the horizontal Laplacian. 

We employ the standard Galerkin procedure to describe the convection patterns at 
the onset of secondary instability in relatively high Prandtl number fluids. The spatial 
dependence of all vertical velocity and temperature field are expanded in a Fourier series, 
which is compatible with the stress-free fiat conducting boundaries and periodic boundary 
conditions in the horizontal plane. Here, we restrict ourselves to standing patterns and, 
hence all time-dependent fourier amplitudes will be assumed to be real. The expansions for 
all the fields are truncated to describe straight cylindrical rolls and patterns arising from 
the nonlinear superposition of mutually perpendicular set of rolls of the same wavenumber. 
Perturbative fields with the same wavelength in mutually perpendicular directions are 
likely to occur in square containers. The vertical velocity and 9 then read as 

V3 — [Wioi{t)cos{kxi) + Woii{t)cos{kx2)]sin{7rx3) + Wu2{t)cos{kxi)cos{kx2)sin{2Trxs) 



d — [Qioiit)cos{kxi) + Qoiiit)cos{kx2)]sin{Trx3) + Qii2{t)cos{kxi)cos{kx2)sin{2Trx3) 



The horizontal components of the velocity field can easily be computed using the equation 
of continuity. Projecting Eq.(4) for V3 and the equation for 9 (Eq. 3) on above modes, we 
arrive at the following minimal mode Lorenz-like model 



+ ..., 



(5) 



+ Qoo2{t)sin{27iX3) + .... 



(6) 






-2ad^S + a^T - -i^X^X2 
d? 2d' 



(9) 



tS 



tT 




(10) 



tZ 



-hZ + X.Y 



(11) 
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where the Unear modes X 



and Y 



^11 t. 2 I '-'101 



\ J^2 J \ '-'oil J 

represent the vertical velocity and the temperature field respectively. The nonlinear 
mode Z = ^^0oo2 denotes heat flux across the fluid layer. Modes S = f^l4^ii2 and 
T = 7%0ii2 are essential to describe nonlinear coupling between two sets of mutually per- 
pendicular rolls. The parameters defined by k'^ = tt^ + A;^, ~ fj' = f?; ^ = ^^gt^^ 
are, in general, wavenumber dependent, while other parameters given hy k'^ — Qc — 



TT^ + A;^, T — \, b — ^ are constants in the model. The reduced Rayleigh number 

Qc Qc 

r = is the control parameter of the problem. We set, hereafter, k — k^-, the wavenum- 

Qc 

ber of stationary rolls at the onset of primary instability. This makes ~ — 1, — ^ 
and 6 = I in our model for convection. 

The steady two dimensional(2D) rolls parallel to Xi-axis is given by Xi = = Yi. 
This makes the nonhnear modes S and T to decouple from the system. Our model then 
reduces to well known Lorenz model with steady solutions given by = ^2 ^6(r — 1) 
and Z = r — 1. However, these solutions become unstable at much lower vahics of reduced 
Rayleigh number r than that predicted by the original Lorenz model[7]. The lower curve 
in the Figure 1 shows the critical value Tq of reduced Rayleigh number, at which the 
2D rolls become unstable, as a function of the Prandtl number a. The critical value To 
increases with increasing a. 



The steady and perfect squares are obtained by setting = X^ in our model. Such 
solutions are given by, 

V2 _v2 \r _ 2o-(r-4d'^)rXi+rXf ^ _ sign{X^X-i) ( g -y^ V _1_ /72 

2ad4X2-2(7(r-4d6)(l+-J-) "^"^ > 

T = ^^g^i|^(aci^X.Y + rXl), Z = i(X.Y), where Xi is a possible real root of the 
algebraic equation 
Xf + AXl + 5 = 

with A = and B = ■'-\^')^r-i) _ ^ot l< r < r„(= Ad^), which defines 

( b +2) ( 6 +2) 

the validity range of the model, B is negative and consequently positive root of Xl exists. 
This implies that the convective patterns in the form of perfect squares exist. However, 
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they are found to be always unstable. 

We study the time-dependent solutions by numerically integrating our model. We 
investigate the dynamics of the patterns as a function of increasing r in two ways- first by 
using the data from previous solution as initial conditions for a new r; and secondly by 
choosing randomly small values of all variables as initial conditions. We get the same result 
with suitably chosen initial conditions. For Vo < r < (see Fig. 1), we find oscillatory 
solutions. The stationary rolls become unstable and a new set of rolls perpendicular of 
the old one develops. The competition between these two sets of rolls leads to a time 
periodic sequence of rolls and patterns arising from nonlinear superposition of two sets 
of rolls. For r > r^^, the solutions become irregular in time indicatiing more modes are 
required to study this regime of parameter space. 

Figure 2 shows the temporal behaviour of the amplitude Xi of the new set of rolls 
as well as the amplitude X2 of the old set of rolls. While Xi oscillates with zero mean, 
X2 oscillates around a finite value. The time period of Xi is always double of that of X2. 
In fact all the modes (Xi, li, S, T) that develop at r > oscillate with same period and 
with zero mean. The periods of I2 and Z are the same as that of X2. The frequency 
/ = l/T of oscillation of X2 at the onset of oscillatory instabihty increases with increasing 
Prandtl number (see Fig. 3). 

The projection of limit cycle on Yi — Y2 plane is shown in Fig. 4. The mean 
of oscillating amplitude of 2D roll mode (1^2) decreases and the perturbative amplitude 
(Yi) increases with increasing r. Figure 5 shows the sequence of shadowgraphs of the 
patterns for < r < rj^. The intensity of old 2D rolls decreases when a new set of rolls 
perpendicular to the old one appears. While the spatial positions for up and down flows 
in the old set remain fixed, these positions alternate for new set of rolls. Consequently, 
the corner positions of square patterns slide back and forth along the set of old rolls. The 
standing waves in form of square patterns do not have four-fold symmetry although they 
preserve inversion symmetry. The time periodic appearances of asymmetric squares, 2D 
rolls, asymmetric squares again with shifted positions of the maximum up or down flows 
by half a wavelength and 2D rolls form a limit cycle. Although the competetion between 
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rolls and symmetric squares is known in binary mixtures[8, 9], but time periodic sequence 
of asymmetric squares and two dimensional rolls in pure fluids is qualitatively new. 

We now test the stability of this limit cycle by introducing the vertical vorticity in 
the model. We expand the vertical vorticity c^s as 

^3 = [Cioi{t)cos{kxi) + (oii(t)cos{kx2)]cos('KX3) + (iio{t)cos{kxi)cos{kx2) 

+ ..., (12) 

Simultaneously, for the consistency of the model, we add a mode W112 in the expansion 
of vertical velocity vs, and a similar mode 6112 in the expansion of temperature field 9. 
This makes a Lorenz-like model consisting of twelve modes. The results of this model for 
fluids with a > 2 exactly reproduces the results of the seven mode model discussed earlier. 
Only for a < 2 and much higher values or r, vertical vorticity is excited. So for a > 2 
and perturbations of wavenumber same as that of 2D rolls, we always see a stable limit 
cycle. The results of our model indicate that for fluids with moderate and large values of 
Prandtl number enclosed preferably in a square container, the time periodic competition 
between the asymmetric squares and rolls should be realizable. 

Acknowledgement: We acknowledge support from DST, India. U. Ghosal, a summer 
trainee from I I T, Kharagpur, acknowledges partial support from I S I, Calcutta. 
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FIGURES 



Figure 1. The region of r and a space showing space rolls, asymmetric squares and 
irregular solutions. The lower curve shows the values of reduced Rayleigh number To 
above 2D rolls become unstable and time periodic asymmetric squares appear. The upper 
curve corresponds to the onset of irregular solutions in the model. 

Figure 2. The convective amplitudes Xi and X2 of new and old sets of rolls respectively 
with respect to time for r = 14.5 and r = 14.7 with value of cr = 10. The time period of 
the new set of rolls is always double that for the old set of rolls, which became stable. 

Figure 3. The plot of the frequency f — 1/T oi the amplitude X2 as a function of 
Prandtl number a at the onset of secondary instability. 

Figure 4. The projection of limit cycle on Yi — Y2 plane for different values of r for 
a = 10. Upper and lower plots correspond to r = 14.5 and r — 14.7 respectively. The 
mean of 2D roll mode (12) decreases and the perturbative amplitude Yi increases with 
increasing r. 

Figure 5. Time-periodic sequence of asymmetric squares and rolls in shadowgraph for 
r = 14.7 and a = 10. Shadowgraphs a,b,c correspond respectively to maximum, zero and 
minimum values of Yi in Fig. 4. 
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